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ABSTRACT 


A fully developed turbulent boundary layer with a zero pressure gradient 
has been studied in a subsonic wind tunnel. The pressure fluctuations on the 
wall of the wind tunnel were measured by small flush-mounted microphones. 
Data for the dimensionless spectral density of the pressure fluctuations, 

P (f)/p2 US 5*, as a function of the frequency parameter f 6*/U, are presented; 
oq Up, 5*, and f are respectively the spectral density, fluid density, free 
stream velocity, boundary layer displacement thickness, and frequency. Data 
for the coefficient (p2)4/\6 p Us are also presented. The transverse cross corre- 
lation for the pressure fluctuations was studied by using two flush mounted 
microphones. Rather than the longitudinal cross correlation, the longitudinal 
cross spectral density was studied, yielding a measure of the coherence of the 
pressure fluctuations at two points as a function of fx/U,, where z is the longi- 
tudinal spacing of the two points and U, is the effective convection velocity of 


the spatial pattern. 


INTRODUCTION 


Pressure fluctuations necessarily coexist with velocity fluctuations in a turbulent 


flow. For incompressible flow the pressure and velocity are related by the equation 


2 
0 u,u,; 


-p [1] . 
0&0 x, 


which can be obtained by combining the continuity and the Navier-Stokes equations. The 
solution to this equation is given in terms of an integration over all space and shows that 
for the case of boundary-layer flow, the vanishing of the velocity fluctuation at the wall does 
not imply the vanishing of the pressure fluctuations at the wall. 

It is these pressure fluctuations on the wall adjacent to the turbulent boundary layer 
that are the focus of interest in this study. Some measurements of the characteristics of 
these pressure fluctuations are reported and the meaning of the measurements are interpreted. 

There are several reasons for being interested in the characteristics of the pressure 
fluctuations. A direct and immediate need exists for data in the field of noise due to unsteady 
flow. When a turbulent boundary layer develops on a flexible wall, the motion of the wall 
which results from the pressure fluctuations generates a sound field. Some of the noise in 


aircraft has been attributed to this mechanism. Kraichnan!’? has given a theoretical 


lReferences are listed on page 12. 


discussion of this problem which begins with the interpretation of Equation [1] for the case 
of boundary-layer flow. 

There are also fundamental reasons for being interested in the pressure fluctuations 
associated with a turbulent boundary layer. There is considerable uncertainty as to how 
energy is produced, convected, diffused, and dissipated in the turbulent boundary layer. In 
particular, the role of pressure fluctuations is but vaguely understood. Although this work 
does not measure the pressure fluctuations inside the boundary layer, it does constitute a 
start on this very challenging problem and it is believed that a further extension of this work 
will prove valuable. 

To obtain the data for this report, the boundary layer on the wall of a subsonic wind 
tunnel was studied. Pressure transducers with a very small active area were mounted flush 
with the wall. Using conventional sound analysis equipment consisting of a narrow tunable 
filter and an rms meter, it was possible to obtain the spectral density” of the pressure fluc- 
tuations at a point on the wall. 

By measuring the pressure fluctuations in the total frequency range, it was possible 
to obtain the mean square value for the pressure fluctuations. This value was also checked 
by integrating over the spectral density. 

By using two pressure transducers and studying the correlation and the cross spectral 
density between the pressure fluctuations at the two points, it was possible to study the 
spatial pattern of the pressure fluctuations and how rapidly it evolves as it is convected 
downstream. By measuring the time delay necessary to maximize the correlation between the 
two points it was possible to measure the velocity at which the pattern was convected. 


These, and other results, are presented and discussed. 


EXPERIMENTAL APPARATUS 
THE WIND TUNNEL 


The wind tunnel was a closed circuit subsonic tunnel. The velocity range of the 
tunnel was from 50 to 200 ft/sec. The transverse dimensions of the working section were 
20 by 15 inches. The measurements were performed on the wall of the working section, 5 feet 
downstream from the entrance nozzle. At this point the displacement thickness of the 
boundary layer was 0.105 inch at a velocity of 100 ft/sec. By comparing the mean velocity 
profiles with those given by Klebanoff? it was judged that the boundary layer was nearly well 
enough developed to show similarity. Unfortunately, the working section was not long enough 
to permit working further downstream so as to prove that the boundary layer was exhibiting 


similarity. 


* 
This term will be defined in the section giving the results of the measurements. 
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Fieure 1 — Schematic of Adaptation of the Altec 21-BR-180 Microphone 
for Spectral Density Measurements 
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Figure 2 — Schematic of Microphone Adapted for Cross Spectral 
and Correlation Measurements 
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THE PRESSURE TRANSDUCER 


The pressure transducer was an adaptation of the Altec 21-BR-180 open-faced 
microphone.* This microphone was modified as shown in Figure 1 to produce the transducer 
used for the spectral density measurements. A different modification was required for the 
correlation measurements. In both cases the upper frequency limit is determined by a 
resonance between the mass of the air in the orifice and the stiffness of the air in the chamber 
adjacent to the diaphragm. Well above this resonance frequency the sensitivity is inversely 
proportional to the frequency range. The resonant response was damped by a cotton plug in 
the orifice so that some of this frequency range was salvaged. 

For the correlation measurements it was necessary to measure at two points separated 
by a small distance. This necessitated the use of the transducer in the form shown in Figure 2. 
Unfortunately, this arrangement resulted in a much lower resonant frequency due to increased 
mass in the tube. Using this transducer, arrangements were made to have a series of holes 
spaced at 0.125 inch apart in the transverse and longitudinal directions which permitted 
correlation measurements to be made up to 1.5 inches in 0.125-inch steps. 

The sensitivity of the transducer was - 61 db below 1 volt for 1 dyne/cm?. For the 


spectral density measurements the usable frequency range was up to 8500 cps and for the 


ow 


correlation measurements the usable frequency range was up to 2000 cps. The lower limit 
on the frequency range was imposed by wind tunnel noise and not the transducer. 

In order to prove that the active area of the transducer was sufficiently small to insure 
that the pressure was coherent over its entire face, various size orifices were used. It was 
found that orifices up to 0.125 inch in diameter were satisfactory for obtaining the data on 
spectral densities. Smaller orifices were used for the correlation measurements in order to 


improve the accuracy in the measurement of the separation distance between the transducers. 


ELECTRICAL EQUIPMENT 


The spectral density of the pressure fluctuations was obtained by the use of the 
Muirhead-Pametrada wave analyzer. A frequency band of a nominal 10 percent was used in 
the analysis. The levels were read on the Ballantine True rms voltmeter. In the section 
giving the results of the measurements, the details of how to obtain the spectral density from 
the observed data are given. 

The correlation between the pressure fluctuations at two points was measured using 


rather simple apparatus. The basis for the method used is contained in the algebraic identity 


€4e> (ec, +e,)? -(; —e5)- 


Sage paaaae 
(Gros 4(e2¢2)% 


where the term on the left is the correlation between ey and Co. The bars indicate time 


averages. e, and e, are here regarded as voltages that are proportioned to the fluctuating 


1 
pressures p, (¢) and p,(¢), respectively, where the subscripts denote the two points at which 
the measurements are performed. With reference to Figure 3, the addition and subtraction of 
the voltages were accomplished by the transformer. The time averages were accomplished by 
the Ballantime True rms voltmeter, whose time constant had been increased to two seconds. 
This produced satisfactory smoothing of the fluctuations. Measuring the correlations is then 
a matter of making four meter readings, [(e, + Ep yal lei Bryce, (2 2)4, and e 2h, and 


doing the computations indicated by the above algebraic identity. 


Wave R.M.S. 
Analyzer Meter 


Figure 3 — Schematic Diagram of the Equipment Used in the 
Measurement of the Cross Spectral Density 
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The wave analyzer shown in Figure 3 serves the function of measuring the correlation 
in a restricted frequency band. As explained later, if observations are made using a suffi- 
ciently narrow frequency band, it is possible to use these measurements to obtain the cross 
spectral density between the pressure fluctuations at two points. 

The convection velocity for the partially frozen pattern of the pressure fluctuations 
was determined by examining the correlation and cross spectral density with time delays. 
These were produced by a delay line consisting of a multistage low-pass filter. This delay 
line gave time delays in 40-microsecond steps up to 6 milliseconds. The delay line had a 


usable frequency range from zero to a little above 5000 cps. 


EXPERIMENTAL RESULTS 
SPECTRAL DENSITY OF THE PRESSURE FLUCTUATIONS 


The spectral density of the pressure fluctuations P(f) is defined as> 
P(f) = ji : |A 2 
= li — 
Fee T rf)! 


where 


T 


An(t)= f p(eniet at 


-T 


is the Fourier transform of the pressure fluctuations over the time 27. 
The basis for the measurement procedure used in obtaining the spectral density begins 


with the fact that the spectral density and the autocorrelation function of p(¢) are Fourier 


transforms. > 


R(t) =p(t)p(t +7) = I P(f) cos wr df 
i) 


It then follows that 


R(0)=p2= | P(f) af 
J 


According to this equation one can experimentally obtain the spectral density by measuring 


the mean square output De of a narrow tunable filter of band with Af and unity gain. 


5 


That is to say, 


This was the procedure used in the work reported here. 


It is desirable to present the measured data in dimensionless form. Since the spectral 


density of the pressure fluctuations has the dimensions of (pressure)? x time, the dimension- 


less spectral density was chosen to be P (f)/p? U3 5* where p, Up, and 6* are the fluid density, 


free steam velocity, and boundary layer displacement thickness , respectively. The 


frequency f was nondimensionalized as f5*/U). 


The dimensionless spectral density curves are shown in Figure 4. The data exhibits 


some lack of similarity over the velocity range studied. As mentioned earlier, the work sec- 


Figure 4 — Spectral Density of the 
Pressure Fluctuations 
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tion of the tunnel was not long enough to 
permit observations further downstream so 
as to prove that the boundary layer was ex- 
hibiting similarity. 


The coefficient 


can be obtained either by integrating over 
the spectral density or by direct measurement 
of D2. In view of the lack of similarity ex- 
hibited in Figure 4, a definite value cannot 
be obtained. The values obtained for K have 
an average value of about 9.5 x 107% which 
in excess of the value obtained by Willmarth.® 


A surprising feature of the spectral 


density is that it does not decrease appreciably at the lowest frequencies measured. For 


these frequencies, which correspond to wave numbers much smaller than the reciprocal of 


the characteristic lengths of the boundary layer, Kraichnan has argued that the spectral density 


should be proportional to the square of the frequency.! There is some evidence that inter- 


mittency of the boundary layer has a strong influence on the pressure fluctuations at the 


lower frequency. Qualitative observations at the lower frequencies revealed an intermittency 


in the amplitude of the spectral density which could not be accounted for by the randomness of 


the process. This intermittency gave pressure versus time records not unlike those observed 


by Klebanoff? for the velocity fluctuations in the boundary layer. This type of behavior 
diminished with increasing frequency and was not noticeable for values of f 5*/U, exceeding 
4x 107%. Further support for the possibility that intermittency plays a role at the lower fre- 
“quencies can be found in the fact that at these frequencies, the measurements of the convec- 
tion velocity yielded a value of about 0.8U). That part of the boundary layer that moves with 
a mean velocity of 0.8U) is strongly intermittent. 

The validity of the above conjectures could be established by repeating this study 
using fully developed pipe flow. 

In view of the preceding, it would have been desirable to have procured data at lower 
frequencies. Unfortunately, this was not possible due to the tunnel noise which limited the 
lowest usable frequencies to between 100 and 200 cps. The signal-to-noise ratio above this 
range was always in excess of 20 db. This figure was obtained by placing the pressure pick- 
up in the middle of the working section of the tunnel where the flow was essentially laminar. 
Actually, some noise was undoubtedly induced by the flow around the microphone, so that 
a better signal-to-noise ratio might actually have existed at the wall. 

Vibration pickup by the transducer was not a problem. By simply plugging up the hole 
in the transducer, the vibration excitation could be measured, which in this case was too low 


to even consider. 


LONGITUDINAL CROSS SPECTRAL DENSITY 


In this work it was decided to study the longitudinal cross spectral density rather than 
the longitudinal correlation. Since this is a novel approach, a few words of explanation are 
needed in order to explain the measurement procedure. 


The cross spectral density of the pressure fluctuations is defined as 
A i 
Pio(f) = ia @ ene e@) 


where the asterisk denotes the complex conjugate, and where 


T 

Ar(th= f pi(tye trae 
=7 
T 

Br(t)= | paltve tat 
i 


with the subscripts 1 and 2 denoting the two observation points. 


Using the arguments of Reference 5 which show that the spectral density and the 
autocorrelation are Fourier transforms, it can be shown that the cross correlation and the 


cross spectral density are also Fourier transforms: 


st iL ee : 
nl) = Pl) PAle son) = il 2 (yer Gh 


—oo 


Since P,,(f) is complex, we can write 
le (i) a Us off) VF WV (fF) 


where it follows from the definition of P,,(f) that the real functions U,,(f) and V,,(f) are 


even and odd, respectively. Accordingly, we can write 


co co 


R,,(7) = ii U,2(f) cos wr df - il ng WSs Ore cey 
0 0 


This last equation provides the basis for an experimental method of measuring the real 
and imaginary parts of the cross spectral density. If the cross correlation is measured with 


zero time delay, then 
RF, (9) i= i U,.(f) af 
) 


If the cross correlation is measured in a narrow band of frequencies Af as shown in Figure 3, 
then 


where p12 (0) denotes the cross correlation measured in the narrow frequency band. The 


real part of the normalized cross spectral density w,,(f) is then 


had) 
La Qe 


Uyo(f) 


Similarly, V,,(f) is obtained by using a time delay 


where f, is here the midband frequency of the narrow filter. The imaginary part of the nor- 


malized cross spectral density v,,(f) is 


Var) 
Oy AGye 


ff) 


It is noted that the measurement of u,,(f) and %,9(f) yields more information than the 
more commonly measured normalized cross correlation since the normalized cross correlation 
may be obtained by integration. 

For this work P, (f) = P, (f) since the two measurement points are sufficiently close 
together that the boundary layer has not changed in character. 

The real and imaginary parts of the cross spectral density as a function of f, x/0.8 U, 
are shown in Figure 5. The distance between the two measurement points has been desig- 
nated as x, and the midband frequency of the narrow filter has been designated as f,. The 
factor of 0.8 will now be accounted for. 

By measuring the cross correlation in a narrow band of frequencies, and by determining 
the time delay T that maximized this correlation, an effective convection velocity U, = x/T 


of the spectral pattern was determined. For all frequencies up to 2000 cps (the upper limit 


Up (f) and vio (f) 
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Figure 5 — Real and Imaginary Parts of the Longitudinal Cross Spectral Density 


of the transducer for this circumstance), it was found that the convection velocity was about 
0.8U,. It was expected that the convection velocity would decrease with increasing values 
of frequency as a consequence of the quadripole nature of the right hand term in Equation [1]. 
The measurement procedure was not sufficiently precise to reveal such a relationship even 
though there is reasonable certainty of its existence. 

It can be seen that the spatial pattern of the pressure fluctuations is not convected 
downstream as though it is frozen. If this were so, the real part of the normalized cross 


spectral density can be shown by a simple computation, for the case of a filter with a rectan- 


‘ ear 
Sa 
fot 0.80, Qnfyz 
= = “= GOS 
“12. \ 0.8U, (2)  \ O8W, 


0.8 U, 


cular frequency response, to be 


where the first factor on the right is an effect due to the finite bandwidth Af. Since, in these 
measurements a constant 10 percent bandwidth was used, this factor is negligibly different 
from unity. 

The departure, for the real part of the data of Figure 5, from cos 27f,z/0.8U, is a 
measure of the development of an uncorrelated component in the pressure fluctuations during 
its travel from the upstream to the downstream measurement point. 

A measure for the uncorrelated component can be formulated mathematically. The 
downstream pressure fluctuation p, (¢) can be resolved into three components whose cross 
spectral densities with the upstream pressure fluctuation are, respectively, real, imaginary, 


and zero. [hat is, 
Py(t) =a(t) + B(d) + y() 


where P, ,=U,>, P} 3 =?Vi25 Pry =0. The vanishing of the spectral density is equivalent 
to saying that the correlation for p, (¢) and y(¢) is zero. This resolution can be affected using 


the definition for cross spectral density. The spectral densities for a(t), 8(¢), and y(¢) are 


2 

ee =U,5/P, 
2 

pelle 


1 


2 / 
Py al Ae, 


10 


Figure 6 — Longitudinal Coherence of the Pressure Fluctuations 


so that 


y2(t) = fejma 


(0) 


and can be computed in terms of the measured data P,(f), P,(f), and |P,,| 2 = Une + Vaceae 


Perhaps of greater interest is the quantity 


2 
(len 
P,P, P, 


which can be taken as a measure of the coherence of Py (t) with P,(¢). This quantity is 
plotted in Figure 6. 

There are several important limitations in the data presented in Figure 5. One trend 
in the data that is desirable to investigate is the effect of the parameter f 5*/U,. Unfortu- 
nately, the transducer used for the cross spectral density measurements had a usable frequen- 
cy range up to only 2000 cps. Consequently, only the flat part of the spectral density, as 
shown in Figure 4, could be investigated. The-upper limit of 2000 cps set an upper limit on 
fd*/U, between 0.083 and 0.33. As can be seen in Figure 4, this does not permit investiga- 
tion of the cross spectral density of the pressure fluctuations in the range where the spectral 
density is rapidly decreasing with frequency. With reference to the velocity fluctuations in 
the boundary layer, the region investigated corresponded to the region where the velocity spec- 
tral density is described by the — 5/3 law where inertia forces dominate. It would be interest- 


ing to procure data for the cross spectral density in the range where viscous forces dominate. 
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THE TRANSVERSE CORRELATION 


The transverse correlation was measured using the equipment of Figure 8 with the 
wave analyzer omitted. The data are presented in Figure 7. Unfortunately, these data have 
not been checked carefully since they were taken on the last day the wind tunnel was avail- 
able for this work. Since we are not concerned here with the convection of a nearly frozen 
pattern, the cross spectral densities were not investigated. For the purpose of a cross check, 
however, it would have been valuable to have obtained these data. Another limitation of the 


data is that its variation of speed was not checked because of limitations on experimental 
time. 
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